Relativistic description of BCS-BEC crossover in nuclear matter 
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Abstract 

We study theoretically the di-neutron spatial correlations and the crossover from superfluidity of neutron Cooper pairs in the 'So 
pairing channel to Bose-Einstein condensation (BEC) of di-neutron pairs for both symmetric and neutron matter in the microscopic 
relativistic pairing theory. We take the bare nucleon-nucleon interaction Bonn-B in the particle-particle channel and the effective 
interaction PK1 of the relativistic mean-field approach in the particle-hole channel. It is found that the spatial structure of neutron 
Cooper pair wave function evolves continuously from BCS-type to BEC-type as density decreases. We see a strong concentration 
of the probability density revealed for the neutron pairs in the fairly small relative distance around 1 .5 fm and the neutron Fermi 
momentum kf n e [0.6, 1.0] fm -1 . However, from the effective chemical potential and the quasiparticle excitation spectrum, there 
is no evidence for the appearance of a true BEC state of neutron pairs at any density. The most BEC-like state may appear at 
kf„ ~ 0.2 firT 1 by examining the density correlation function. From the coherence length and the probability distribution of 
neutron Cooper pairs as well as the ratio between the neutron pairing gap and the kinetic energy at the Fermi surface, some features 
of the BCS-BEC crossover are seen in the density regions, 0.05 far 1 < kf„ < 0.7 (0.75) fin , for the symmetric nuclear (pure 
neutron) matter. 
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Pairing correlations, as collective phenomena found in vari- 
ous systems as liquid 3 He, superconductors, atomic nuclei, and 
ultracold atomic gases, play a crucial role in the fermion sys- 
tems. When an attractive interaction between two fermions is 
weak, i.e., the pairing gap is much smaller than the Fermi en- 
ergy, partner fermions can be treated as a delocalized Cooper 
pair with large overlap as described in the Bardeen-Cooper- 
Schrieffer (BCS) theory [1], displaying a strong correlation in 
the momentum space. If an interaction is sufficiently strong, 
fermion pairs become a spatially compact bosonic bound state 
and undergo Bose-Einstein condensation (BEC), showing a strong 
correlation in the coordinate space. As the strength of attractive 
interaction between two fermions increases, the pairing phe- 
nomenon evolves from BCS to BEC state. The transition be- 
tween them, often referred as BCS-BEC crossover, brings a 
new insight into the pairing phenomenon. Although the BCS 
and BEC limits are physically quite different, the evolution be- 
tween them was found to be smooth and continuous J2, S 0] • 

In nuclear physics, neutron pairs are expected to be strongly 
correlated at some density. There are rather general reasons to 



'Corresponding author 

Email addresses: tokiOrcnp . osaka-u .ac.jp (Hiroshi Toki), 
meng j Spku .edu.cn (Jie Meng) 



expect that di-neutron correlations in low-density nuclear mat- 
ter should be significant. It is well known that the bare nucleon- 
nucleon interaction in the 'So channel leads to a virtual state 
around zero energy characterized by the large negative scatter- 
ing length a w -18.5 ±0.4 fm [5]. This implies a very strong at- 
traction between two neutrons in the spin singlet state, although 
the neutron-neutron interaction is not strong enough to form a 
two-body bound state in free space. Furthermore, theoretical 
predictions suggest that at around 1/10 of the nuclear satura- 
tion density po, the 'So pairing gap may take a considerably 
larger value than that around po Ma [ZD . 

Weakly bound neutron-rich nuclei provide optimum circum- 
stance to study di-neutron correlations in the low-density re- 
gion, where the density around surface is unsaturated and the 
c oupling s to the continuum spectra play an essential role JH 
SOU HI 12]. Originally the possible existence of a di-neutron 
near the surface of nuclei was predicted by Migdalin 1972 [13]. 
In his work, it was shown that although di-neutron is unbound in 
vacuum, it may form a bound state in a potential well if there is 
a single-particle level with energy close to zero. Subsequently, 
plenty of theoretical investigation based on either schematic or 
microscopic models displayed a clear signature of the strong di- 
neutron correlation in the ground state, not only in l igh t closed- 
shell core plus In or An nuclei lfl4l[8T[lHfl6l[H[l2l[l8i[l9ll, but 
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also in medium or heavy neutron-rich nuclei 1 2(J 21 , 22, 231. as 
well as in infinite nuclear matter [24, 25]. It has been revealed 
that the concentration of small sized Cooper pairs in the surface 
of superfiuid nuclei is a quite generic feature 1I23I1 . In fact, such 
a di-neutron correlation is in connection with the BEC-like be- 
havior in infinite nuclear matter at low densities. Recently, a 
strong low-lying dipole strength distribution of the two-neutron 
halo nucleus 11 Li has been observed experimentally [26]. The 
spectrum was reproduced well bya three-body model with a 
strong di-neutron correlation 11271 12811 . 

Theoretical and experimental progress on di-neutron corre- 
lations in weakly bound nuclei has stimulated lots of renewed 
interests in possible BCS-BEC crossover of di-neutron pairs. 
The systematical analysis was performed in nuclear matter by 
using the bare force given by a superposition of three Gaus- 
sian functions, the finite-range Gogny interaction, and the zero- 
range contact interaction without or with medium polarization 
effects i29l, [30I 31]. It was shown that the correlations be- 



tween two neutrons get large as density decreases, and the spa- 
tial structure of neutron Cooper pairs changes. As a result, the 
BCS-BEC crossover occurs in the uniform matter at low den- 
sities. In particular, for the screened pairing interaction, a di- 
neutron BEC state is formed in symmetric matter at around 
0.00 1 pd. The above results have been confirmed by a further 
study, where the change of sign of the density correlation func- 
tion at low momentum transfer is taken as a criterion of the 
BCS-BEC crossover rf32ll . In addition, the coexistence of BCS- 
and BEC-like spatial structures of neutron pairs was studied in 
the halo nucleus 11 Li as well 03311 . It has been shown that as 
the distance between the center of mass of the two neutrons and 
the core increases, the two-particle wave function changes from 
the weak coupling BCS regime to the strongly correlated BEC 
regime. This result clearly confirms that a strong di-neutron 
correlation between the valence neutrons is present on the sur- 
face of the nucleus. 

Details of the effective nucleon-nucleon force in nuclei are 
as yet not completely clarified since most of the experimental 
data in nuclear spectroscopy are not very sensitive to the details 
of the nucleon-nucleon force. Hence, in most of investigations 
on nuclear pairing correlations, for convenience, the effective 
nuclear forces such as zero-range contact forces or finite-range 
Gogny force, are used in the particle-particle (pp) channel. The 
effective interactions in relativistic mean-field (RMF) theory 
are also used in the pp channel 134, 35]. In order to obtain 
reasonable values for the gap parameter, one has to introduce 
an effective factor [35J. Hence, the properties of Cooper pairs 
and BCS-BEC crossover from these existing phenomenologi- 
cal calculations need further check. Furthermore, in the low- 
density limit, the interaction in the pp channel approaches the 
bare nucleon-nucleon interaction. Therefore, it is interesting 
to use the realistic bare nucleon-nucleon interaction in the pp 
channel and explore the BCS-BEC crossover at low densities. 

In this work, the di-neutron spatial correlations in the 1 Sq 
channel will be studied in a relativistic pairing theory, i.e., rela- 
tivistic Hartree-Fock-Bogoliubov (RHFB) theory [34], with the 
realistic bare nucleon-nucleon interaction, i.e., the relativistic 
Bonn potential l36Tl . in the pp channel. Then the BCS-BEC 



crossover phenomenon and the possibility of di-neutron BEC 
states at low densities will be discussed. 

The RMF theory ll3~7ll has attracted lots of attentions dur- 
ing the last decades due to its great success in describing both 
nuclear matter and finite nuclei near or far from the stability 
line 0|i]. In the RHFB theory, meson fields are treated 
dynamically beyond the mean-field theory to provide the pair- 
ing field via the anomalous Green's functions 113411 . In the case 
of infinite nuclear matter, the Dirac-Hartree-Fock-Bogoliubov 
equation reduces to the usual BCS equation. For the l So chan- 
nel, the pairing gap function A{p) is, 



v np( k 'P)-^r-k 2 dk, 



(1) 



where v pp (k, p) is the matrix elements of nucleon-nucleon in- 
teraction in the momentum space for the So pairing channel, 
and Ek is the quasi-particle energy, 



E k = ^(s k -p) 2 +A(k) 2 
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with the single-particle energy e*, and the chemical potential p. 
The corresponding normal and anomalous density distribution 
function have the form, 
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2E k ' 



(3) 



The single-particle energy s k follows from the standard RMF 
approach 03711 . 



En + Vfc 2 + M* 2 , 



(4) 



where the scalar mass M* = M + £s , So the vector potential and 
S5 the scalar potential. For nuclear matter with given baryonic 
density pi, and isospin asymmetry £ = (p n - p p )/pb, the above 
equations can be solved by a self-consistent iteration method 
with no-sea approximation. 

The relativistic Bonn potential is used in the pp channel, 
which has a proper momentum behavior determined by the scat- 
tering data up to high energies IBal . It is defined as the sum of 
one-boson-exchange (OBE) of the several mesons <p = <t,oj, tc, p,rj,6. 
The matrix elements v pp (k,p) is 



v pp (k,p) = V -H±-A <p (k,p)D <p (q 2 )Fl(q 2 ), 



where e* k is the effective single-particle energy 



e* = V/fc 2 + M* 2 . 



(5) 



(6) 



D^iq 2 ) is a meson propagator with momentum transfer q = k- 
p, and F,f,(q 2 ) is a form factor of monopole type, 



1 



q 2 +m 2 q 2 + A 2 



(7) 



with the mass of meson, and the cutoff parameter. rj^, and 
A<f,(k,p) are the vertex functions of the relativistic Bonn poten- 
tial. For the 'So pairing channel, the matrix element v pp (k, p) is 



2 



just the average of v pp (k,p) over the angle 6 between the vec- 
tors k and p, 



v PP (k, p) 



r v p P 

Jo 



(k, p) sin 6dG. 



(8) 



Bonn-B potential [36] is adopted for v pp (k,p). For the mean- 
field calculation in the particle-hole (ph) channel, the effective 
interaction PK1 j40ll is used, since the results do not depend 
sensitively on various other parameter sets [41]. 
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Figure 1: Neutron Cooper pair wave function in coordinate space 'P p!1 j 1 (r) as a 
function of the relative distance r between the pair partners at several neutron 
Fermi momenta kf „ in symmetric nuclear matter = 0, black solid lines) and 
pure neutron matter (f = 1, red dashed lines). Notice the different scales on the 
ordinate in the various panels. 

To investigate the spatial structure of neutron Cooper pairs, 
it is useful to look into its wave function represented as a func- 
tion of the relative distance r = \r\ — r2\ between the pair part- 
ners. The Cooper pair wave function in momentum space *F pa ii (k) 
is just the anomalous density k^, and its coordinate representa- 
tion is deduced from the Fourier transform, 



YpairW 



(2tt) 3 / 



K k e 



r dk, 



(9) 



where C is the constant determined from the normalization con- 
dition 



2 r 2 dr-. 



1. 



(10) 



In Fig. \T\ the wave function of the neutron Cooper pairs 
^paiiM is shown as a function of the relative distance r be- 
tween the pair partners at different neutron Fermi momenta kf n . 
Nearly identical results are given for the symmetric nuclear 



matter = 0) and pure neutron matter = 1), except for 
kf„ = 1.2 fm -1 , where a larger amplitude of the first peak is 
obtained in pure neutron matter. 

The radial shape of ^Ppaiif/) changes as density decreases. 
When kf„ = 1.2 fm -1 and 0.8 fm -1 , the profile of the wave 
function T^ak (/) shows an oscillatory behavior convoluted by a 
decreasing exponent, which is a typical behavior of BCS state. 
A significant amplitude outside the first node is observed. As 
density goes down to kf n = 0.4 fm -1 and 0.2 fm -1 , the wave 
function becomes compact in shape and the oscillation disap- 
pears, resembling the strong coupling BEC-like bound state. 
This indicates that a possible BCS-BEC crossover may occur 
in uniform matter at such low densities. At very dilute density 
of k F „ = 0.05 fm -1 , the wave function becomes more extended 
again and approaches to zero slowly. 

Turning ones attention to short distance r, a suppressed am- 
plitude around r = fm is displayed at all densities, that is at- 
tributed to the strong repulsive contribution of v pp (k,p) at high 
momentum. In addition, for all the densities a peak appears 
around r = 1.0 fm, around which it has been confirmed that 
the pairing potential in the coordinate space v pp (r) reaches the 
strongest attraction 1 2511 . 

Furthermore, it is noteworthy that the evolution of the wave 
function with the density in Fig^[T]is similar to those for the 
proton-neutron pairing in ref. [42] in which a true bound state 
of deuterons is predicted in the low-density limit. Thus, one 
needs to explore whether a true di-neutron BEC bound state 
appears in nuclear matter as well. 
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Figure 2: A two-dimensional plot for the probability density r 2 |f ' pa i r (r)| 2 of the 
neutron Cooper pairs as a function of the neutron Fermi momentum kp „ and the 
relative distance r between the pair partners in symmetric nuclear matter. 

To illustrate a clear picture of di-neutron correlations, a two- 
dimensional plot for the probability density ^I^Ppa^r)! 2 of the 
neutron Cooper pairs for symmetric nuclear matter is shown 
in Fig. 12 as a function of the neutron Fermi momentum kp n 
and the relative distance r between the pair partners. A strong 
concentration of the probability density is revealed for the pair 
partners in a surprisingly small relative distance of r ~ 1 .5 fm 
in the density region about 0.6 firT 1 < kf n < 1.0 firT 1 . This 
behavior has the same physical essence as the di-neutron con- 
figuration in the surface of superfluid nuclei as discussed in 
Refs. J22, 23, 33]. As a matter of fact, the pairing gap at the 
Fermi surface Af„ = A(kf„) also depends strongly on kf„, and 



3 



achieves a large value together with probability density in such 
a density region. 

When the Fermi momentum is close to the saturation den- 
sity, the probability density ^PFpairfV)! 2 indicates a considerable 
distribution outside the first peak and extends to a large distance 
r. This phenomenon is well understood by the orthogonaliza- 
tion of the wave function for the paired particles with the re- 
maining particles due to the Pauli principle, which expels the 
Cooper pair wave function outside. As density decreases, these 
amplitudes die out gradually and a compact structure appears at 
about 0.2 firT 1 < kf„ < 0.5 far 1 . While going to lower densi- 
ties of kf n < 0.2 firT 1 , the probability density has a very small 
value and exhibits a long tail at large relative distance r. This 
behavior may be responsible for the appearance of halo struc- 
ture in several neutron-rich nuclei at far distance away from the 
center UEHi. 

Since a strong di-neutron correlation is revealed at low den- 
sities from the Cooper pair wave function and the correspond- 
ing probability density, it is deserved to study the possibility of 
Bose-Einstein condensation and the BCS-BEC crossover phe- 
nomenon of di-neutrons as density decreases. 
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Figure 3: Effective neutron chemical potential v„ as a function of the neutron 
Fermi momentum kf „ in symmetric nuclear matter (black solid line) and pure 
neutron matter (red dashed line). 

In the whole debate around BCS-BEC transition the chemi- 
cal potential p plays a central role. It has been proved in an easy 
way that the gap in Eq. (Q]) together with the normal and anoma- 
lous density distribution function in Eq. (0 goes over into the 
Schrodinger-like equation for the neutron pair wave function 
Tpaj^fc), and the corresponding energy eigenvalue is related to 
the chemical potential yfl. In the case of the RMF approach, it 
is expressed as 

2e(k)% dir (k) + (1 ~^ Pk) j v pp {k,p)p 2 dpV vm (p) 

= 2v„Y pair (&), (11) 

where e(k) is the neutron kinetic energy, 

e(k)= y/k 2 + M* 2 -M*. (12) 
v„ is the effective neutron chemical potential obtained by de- 



ducting momentum independent part from the chemical poten- 
tial p„, 

v n =p n -L -M\ (13) 

which could be regarded as half of the "binding energy" of a 
Cooper pair. At transition from the BCS regime of large over- 
lapping neutron Cooper pairs to the BEC regime with true di- 
neutron bound states, v„ is supposed to turn from positive to 
negative values. 

In Fig. |3] the effective neutron chemical potential v„ as a 
function of the Fermi momentum kf„ in nuclear matter is given. 
The values of v„ decrease monotonically in both symmetric nu- 
clear matter and pure neutron matter as density goes down, and 
come very close to zero at dilute area but never turn negative. It 
is then expected in this case that a true di-neutron bound state 
does not occur in nuclear matter, but the transition to BEC is 
very close in the low-density limit. 
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Figure 4: Neutron quasiparticle excitation energy Ek as a function of the ratio 
of the neutron momentum to the Fermi momentum k/kp,, at different neutron 
Fermi momenta kf„ in symmetric nuclear matter. 

The qualitative changes in the quasiparticle excitation spec- 
trum have been suggested as another clear distinction between 
the BCS and BEC limits [43], i.e., the momentum correspond- 
ing to the minimum in the excitation spectrum will shift from 
finite momentum in the BCS limit to zero momentum in the 
BEC limit, and the value of the pairing gap also changes from 
A in the BCS limit to yjp 1 + A 2 in the BEC limit. 

In Fig.|U the quasiparticle excitation spectra Ek is shown as 
functions of k/kf„ at several neutron Fermi momenta in sym- 
metric nuclear matter. The minimum of Ek always appears 
around the Fermi momentum kf„ for all the densities, i.e., the 
minimum of Ek approaches to zero momentum with the corre- 
sponding Fermi momentum but become zero only when kf„ = 
0. Similar results are obtained for the pure neutron matter as 
well. Therefore, from the quasiparticle excitation spectrum, 
there is no evidence for the appearance of a true BEC bound 
state of neutron pairing with the realistic bare nucleon-nucleon 
interaction Bonn-B at any density, in agreement with that from 
the effective chemical potential v„. 

The analogies with other low-density Fermi systems, such 
as ultracold atomic gases, could be used to study the BCS-BEC 
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Figure 5: Zero-momentum transfer density correlation function D(q = 0) and 
its normal and anomalous density contributions, 7 p (0) and 7 r (0), as a function 
of the neutron Fermi momentum kp„ in symmetric nuclear matter (black solid 
lines) and pure neutron matter (red dashed lines). 



crossover in nuclear matter as well. The neutron density corre- 
lation function is such a quantity that gives prominence to the 
relative strength between the mean field and the pairing field, 
which is defined as 144 13211 . 



D(q) = I K {q) ~ h{q)- 



(14) 



At zero momentum transfer, the normal and anomalous density 
contributions I p and I K respectively read 



i r°° 
= — pl k2dk > 

^Pn Jo 
= f l4k 2 dk, 

n 2 p n Jo 



4.(0) 
4(0) 

and they satisfy the sum rule 



4(0) + 4(0) = 1. 



(15) 
(16) 

(17) 



The change of sign of the density correlation function at 
low momentum transfer is taken as a criterion of the BCS-BEC 
crossover [44], although further study argued that this criterion 
can be trusted only at small isospin asymmetry 132 1. 

In Fig. [5] the density correlation functions at zero-momentum 
transfer D{0) and its normal and anomalous density contribu- 
tions, I p (0) and I K {0), are shown as functions of kf„ in nuclear 
matter. With decreasing density, I p (Q) first drops down from 
1.0, reaches a minimum at kf„ ~ 0.2 fm -1 (p„/po ~ 10 -3 ), 
and then increases to 1.0 again. While 4(0) gives an oppo- 
site trend and has no contribution when approaching to either 
the saturation density or the dilute density. This leads to a 
peak for the density correlation function D(Q) around kf„ = 
0.2 fm -1 , where the most BEC-like state may appear. It is con- 
sistent with the result of v t , pa i r (r) shown in Fig. Q] where the 
wave function behaves like a spatially compact bound state at 
kf„ = 0.2 fm -1 . Because the anomalous density contribution 



4(0) is always smaller than the normal one / p (0), the density 
correlation function never changes sign at all. According to the 
criterion mentioned above, di-neutrons are not in BEC state but 
just in the transition region between BCS and BEC regimes, 
that verifies the conclusion from the effective chemical poten- 
tial and the quasiparticle excitation spectrum. 
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Figure 6: Probability P(r) for the partner neutrons correlated within r = 3 fm 
(upper-left panel) and the average inter-neutron distance r = d„ (lower-left 
panel) as a function of the neutron Fermi momentum kf„. The black line with 
filled circle is for the symmetric nuclear matter and the red line with filled 
square is for the pure neutron matter. The corresponding rms radius f lms of the 
neutron Cooper pair is plotted in upper-right panel in comparison with the av- 
erage inter-neutron distance d„ (the dashed curve). The ratio A/?„/e/r„ between 
the neutron pairing gap at the Fermi surface Ap„ and the neutron Fermi kinetic 
energy ep n as a function of kf n is plotted in lower-right panel. The referred 
BCS-BEC crossover regions (yellow grid) and the unitary limit (dashed line) 
from the regularized pairing model [29, 30] for P(d„) and A/r„/e/r„ are given 
respectively. 

To describe the BCS-BEC crossover quantitatively, lots of 
characteristic quantities have been introduced [29, 30]. The rms 
radius of the neutron Cooper pair £ rms , i.e., coherence length, is 
a straightforward measure to evaluate the pairing size, 



(r 2 ) 



(18) 



which can be calculated either from the Cooper pair wave func- 
tion ^pairW or from in the momentum space. The latter one 
is more suitable since the gap equation is solved in momentum 
space. The coherence length £ rms has a solid meaning even in 
the strong coupling BEC case or the crossover region between 
BCS and BEC. For comparison, the average inter-neutron dis- 

— 1/3 

tance d n =p n is taken into account. If £ rras > d„, the neutron 
Cooper pair is interpreted as an extended BCS-like pair, other- 
wise it will be considered as a compact BEC-like pair. 

The coherence lengths £ lms of the neutron Cooper pairs as 
a function of the neutron Fermi momentum kf„ in both sym- 
metric nuclear matter and pure neutron matter are drawn in 
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the upper-right panel of Fig. [6] While kp n decreases, the co- 
herence lengths shrinks to ~ 5 fm at density of 0.4 far 1 < 
k F „ < 0.8 far 1 (0.4 far 1 < k F „ < 0.9 far 1 ) for symmetric 
nuclear (pure neutron) matter. Then it expands again at very 
low densities. These features can be understood as follows. 
In comparison with the virtual S -state in the vacuum, the neu- 
tron pairs feel an extra binding due to the pairing correlation 
as density increases; while at large kf„, the neutron pairs will 
strongly increase in size due to the Pauli principle. The dif- 
ference of the coherence lengths between ( — and ( = 1 at 
high densities is caused by the difference of the scalar mass M* . 
When kf„ < 0.7 far 1 (kf„ < 0.75 far 1 ) for symmetric nuclear 
(pure neutron) matter, it is seen that the coherence length £ rms is 
smaller than the average inter-neutron distance d„, which could 
be treated as the signature of the BCS-BEC crossover. 

In order to characterize the spatial correlation of the neutron 
Cooper pairs, one can introduce the probability P{r), 

P(r)= f \V^{r')\ 2 r' 2 dr' , (19) 
Jo 

which gives the probability for the pair partners within a relative 
distance r. 

The probability P(r) for the neutron pair partners correlated 
within r = 3 fm is shown in the upper-left panel of Fig. [6] as a 
function of kp„ in symmetric nuclear matter and pure neutron 
matter. It reaches the maximum value of about 0.75 around 
kp n = 0.75 fm -1 , and goes down at both lower or larger Fermi 
momentum. In the lower-left panel of Fig. [6] the probability 
P(d„), i.e., P(r) for the pair partners within the average inter- 
neutron distance r = d n , is plotted as well. The values of P(d„) 
increase monotonically with decreasing density, and approach 
to 1.0 at dilute area. From a regularized pairing model l45ll . 
it is suggested that the BCS-BEC crossover is determined by 
P(d„) > 0.807 1291 B On . This limits the Fermi momentum kf„ < 
0.7 far 1 in symmetric nuclear matter and kp„ < 0.75 far 1 in 
pure neutron matter for the realization of the BCS-BEC crossover. 

The BCS-BEC crossoveris also explored by the ratio A/r„/e F „ 
between the neutron pairing gap at the Fermi surface Ap„ and 
the neutron Fermi kinetic energy ep n = e(kf n ). If the ratio is 
large enough, the neutron pairing is expected to be in the strong 
coupling regime. From the regularized model, the BCS-BEC 
crossover region is estimated as 0.21 < A n lep n < 1.33 and the 
unitary limit is given at A„/ep„ = 0.69 lE^BOll . 

In the lower-right panel of Fig. [6] the ratios A Fn je Fn plot- 
ted as a function of the neutron Fermi momentum are shown 
for both symmetric nuclear matter and pure neutron matter. As 
density decreases, the values grow up first and go into BCS- 
BEC crossover region at kp„ ~ 0.75 far 1 for symmetric nuclear 
matter and kp„ ~ 0.8 far 1 for pure neutron matter. It is found 
that the curves never intersect with the line of the unitary limit, 
namely, the neutron pairing does not go into the BEC regime. 
Then at very low densities with kp„ < 0.05 far 1 , the results 
decrease again and return back to the BCS region. Based upon 
the above analyses, the BCS-BEC crossover is marked in the 
density region with 0.05 far 1 < kp„ < 0.7 far 1 for the sym- 
metric nuclear matter and 0.05 far 1 < kp„ < 0.75 far 1 for the 
pure neutron matter. 



In conclusion, the di-neutron spatial correlations and the 
BCS-BEC crossover phenomenon for nuclear matter in the l So 
channel has been investigated based on the microscopic rela- 
tivistic pairing theory with the effective RMF interaction PK1 
in the ph channel and the realistic bare nucleon-nucleon interac- 
tion Bonn-B in the pp channel. It is found that the spatial struc- 
ture of neutron Cooper pair wave function evolves continuously 
from BCS-type to BEC-type as density decreases, and a strong 
concentration of the probability density is revealed for the pair 
partners in the fairly small relative distance around 1 .5 fm and 
the neutron Fermi momentum kp„ e [0.6, 1.0] far 1 , which is 
correlated with the di-neutron configuration in the surface of 
superfluid nuclei. In light of the evidence from the effective 
chemical potential, the quasiparticle excitation spectrum and 
the density correlation function, a true di-neutron BEC bound 
state does not occur at any density in nuclear matter. Neutron 
pairing is just in the transition region from BCS to BEC regime 
at low densities, and the most BEC-like state may appear at 
kp„ ~ 0.2 far 1 . From several characteristic quantities, such as 
the coherence length £ rms , the probability P(r) with r — 3 fm 
and r - d„, and the ratio Ap n /e Fn , the BCS-BEC crossover is 
found in the density region about 0.05 far 1 < k Fn < 0.7 far 1 
for symmetric nuclear matter and 0.05 far 1 < kp„ < 0.75 far 1 
for pure neutron matter. The results reveal a strong correlation 
of neutron pairs at such densities. 
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